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ABSTRACT 


An  expression  for  the  kinetic  energy  of  an  in¬ 
finite  liquid  mass  due  to  small  motions  of  a  sub¬ 
merged  cylinder  is  obtained  analytically.  VRien  the 
motion  is  described  in  terms  of  generalized  coordi¬ 
nates  ^  the  corresponding  inertia  coefficients  appear 
in  the  kinetic-energy  expression.  The  analysis  is 
carried  out  for  the  problem  of  a  floating  flexible 
rectangular  box.  Calculations  are  performed  for 
generalized  coordinates  representing  rigid-body  mo¬ 
tion  and  plastic  deformation. 


I. 


Introduction 


The  purpose  of  this  paper  is  to  provide  the  background 
for  the  treatment  of  transient  problems  of  elastic  or  plastic 
box  structures  resembling  a  surface  ship  and  floating  on  the 
surface  of  a  semi- infinite  fluid.  Whenever  the  transient 
problem  is  such  that  the  fluid  may  be  considered  incompres- 
sible>  the  problem  can  be  formulated  in  terms  of  generalized 
coordinates;  the  presence  of  the  fluid  then  produces  only 
terms  which  can  be  derived  from  the  kinetic  energy  of  the 
fluid,  expressed  in  the  generalized  velocities.  The  results 
of  this  paper  are  intended  to  be  applied  to  the  determination 
of  the  elastic  or  plastic  response  of  a  ship-like  structure 
to  initial  velocities  imparted  to  it  by  an  underwater  explo¬ 
sion.  Appropriate  velocity  distributions  may  be  obtained  by 
experimental  or  analytical  means;  an  analytical  approach  has 
been  presented  in  Ref.  1. 

The  present  problem  becomes  amenable  to  analysis  if  it 
is  treated  as  a  two-dimensional  one,  that  is,  if  the  box  is 
regarded  as  an  infinitely  long,  flexible  rectangular  cylinder 
vdiose  boundary  undergoes  small  motions  perpendicular  to  the 
generatrices;  the  latter  remain  straight  and  parallel  to  the 
axis.  The  fluid  is  assumed  at  rest  at  large  distances  frcxn 
the  cylinder,  and  the  motions  of  the  cylindrical  boundary 
are  small  deviations  from  its  rest  configuration. 
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The  problem  may  thus  be  viewed,  approximately,  as  a  bound¬ 
ary-value  problem  in  plane  ideal- fluid  flow  with  prescribed 
normal  velocities  along  a  stationary  boundary.  It  can  there¬ 
fore  be  attacked  by  the  classical  methods  of  plane  potential 
theory,  and,  in  particular,  by  the  conformal  mapping  of  the 
given  boundary  curve  into  a  circle,  and  subsequent  contour  in¬ 
tegration.  The  analysis  is  carried  out  for  the  general  case 
in  Section  II,  resulting  in  the  general  expression  for  the 
inertia  coefficients.  The  expression  takes  the  form  of  an 
infinite  series  whose  terms  depend  on  the  Fourier  coefficients 
of  the  transformed  normal-velocity  distribution.  These  Fourier 
coefficients  may  be  obtained  by  quadrature  (nuxnerical,  in 
general),  provided  the  mapping  function  is  known.  The  mapping 
function  appropriate  to  the  present  case  of  a  floating  rectan¬ 
gular  box  is  developed  analytically  in  Section  III,  and  is 
tabulated  at  the  end  of  the  report  in  a  form  suitable  for  in¬ 
tegration,  for  the  two  beam-draft  ratios  of  2  and  3. 

Calculations  were  actually  carried  out  for  the  special 
case  of  motion  in  which  the  degrees  of  freedom  considered  were 
(i)  the  three  rigid-body  motions  (horizontal  and  vertical 
translations,  and  rotation),  and  (ii)  three  deformations  of 
the  box  representing  possible  plastic  damage.  These  degrees 
of  freedom  are  illustrated  in  Figure  1.  The  results  are  tabu¬ 
lated  in  the  form  of  inertia  matrices  for  the  two  beam-draft 
ratios  considered. 
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An  idea  of  the  accuracy  of  the  calculations  (which,  be¬ 
sides  numerical  integrations,  involve  truncations  of  infinite 
series)  may  be  had  by  virtue  of  the  fact  that  for  the  two 
translatory  degrees  of  freedom  the  inertia  coefficients  may 
be  calculated  exactly.  This  is  shown  in  the  Appendix. 

The  results  of  this  paper  are  also  applicable  to  vibra¬ 
tion  problems,  provided  the  frequencies  are  low  enough  to 
permit  the  assumption  of  incompressibility. 

II .  Virtual  Inertia  Terms  for  a  Closed  Cylinder 

The  kinetic  energy  of  an  incompressible  fluid  in  irrota- 
tional  two-dimensional  flow  in  the  x  ^  plane  is,  per  unit 
distance  in  the  z  direction, 

T  .  I  p  ^  0  d^  ,  (1) 

Where  p  is  the  (constant)  fluid  density,  0  the  velocity 
potential,  ^  the  stream  function,  and  the  integration  is 
carried  out  in  the  positive  direction  aro\and  the  boundary. 
(See,  for  example.  Ref.  2,  p.  66.)  If  the  fluid  is  bounded 
internally  by  a  closed  cylinder,  then,  in  accordance  with  the 
usual  convention,  the  direction  of  integration  is  clockwise 
around  the  projection  C  of  the  cylinder.  On  introducing  the 
complex  potential  function, 

w  «  0  +  i^  , 
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one  may  write  Equation  (l)  as 

T  =  -^  p^w*dw  ,  (2) 

(the  asterisk  denotes  the  complex  conjugate),  since 

w*dw  =  {<t>  -  ii/)  {d<ti  +idil/)  =  2i0d^  + -^  d(0®  +  -  2i0^) , 

and  the  integral  of  a  perfect  differential  around  a  closed 
curve  is  zero. 

The  X  and  ^  components  of  fluid  velocity,  u  and  v, 
are  given  by  the  relation 

u  -  iv  =  dw/dz  ,  (3) 

where 

z  »  X  +  iy  . 


The  normal  velocity  (considered  positive  into  the  fluid)  at 
a  point  z  »  z^  on  the  curve  C  is  therefore 


u 


V 


dy  dx 
"  di  -  dS 


«  lin(dw/ds) 


z>z 


c 


(4) 


where  ^  I  dz  |)  is  an  element  of  arc  length  on  C,  meas¬ 
ured  positive  counterclockwise. 

If  the  z  plane  is  now  mapped  conformally  into  a  C 
plane  such  that  the  outside  of  C  is  mapped  into  the  outside 
of  the  unit  circle. 
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then  for  a  fluid  which  is  at  rest  at  infinity  the  complex 
potential  function  may  be  expressed,  in  general,  in  the  form 


w 


00 

T. 

n=l 


a  +  ib 
n  n 

n 


(5) 


where  the  coefficients  a„  and  b^  are  real.  The  normal 

n  n 

velocity  on  C  is  now  given  by 


u 


V 


d0 

ds 


Z 

n*l 


(a  cos  n0  +  b^  sin  nG)  . 
'  n  n 


(6) 


Consequently,  by  Fourier's  theorem. 


Tlie  arc  length  on  C,  measured  countercloclcwise  from  the 
point  corresponding  to  0  «  0,  may  be  treated  as  a  function 
of  0,  defined  by 

•(9)  -  /  ^ 

vdiere 

(9) 

exp(i0) 

In  particular, 

8(2Tr)  2  p  , 
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the  perimeter  of  C.  We  may  further  define  a  dimensionless 
variable j 

a  =  b/p  , 


and  consider  B  a  function  of  a.  Equation  (7)  may  then  be 
written  alternatively  as 


cos  ne(a) ^ 
sin  n0(o)  J  * 


(10) 


When  the  expression  (5)  for  the  complex  potential  func¬ 
tion  is  used  in  Equation  (2)  and  the  integration  is  carried 
out  along  ^  »  exp(id);  the  resulting  expression  for  the  ki¬ 
netic  energy  is 

T  »  i  ttp  L  --n  .  (11) 

^  n«l  " 

A  small  normal  displacement  of  C  may  in  general  be  re¬ 
garded  as  the  superposition  of  normal  displacements  correspond¬ 
ing  to  different  modes  of  motion  (or  degrees  of  freedom) ,  each 
of  which  may  be  expressed,  at  any  time  t,  as 

fl(o)q^(t)  , 

where  f^{o)  describes  the  shape,  while  q^Ct)  gives  the  am¬ 
plitude.  (it  should  be  noted  that  if  q^  is  a  length,  f^ 
is  dimensionless,  while  if  q^^^  is  an  angle,  f^  has  the  di¬ 
mension  of  length.)  The  general  normal  velocity  is  consequently 
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(12) 


Uy  =  ^  fi(o)qi  . 

The  Fourier  coefficients  a  ,  b  ,  are  then 

n  n 


where 


da 


(13) 


(14) 


The  amplitudes  may  now  be  treated  as  generalized  coordi¬ 

nates  of  the  motion,  so  that  the  kinetic  energy  becomes  a 
quadratic  form  in  the  generalized  velocities: 

T  »  I  2  2  niij  q^q.  ,  (15) 

where 


o  00  c .  c  •  +  a .  S' 

„  .i!Ef  2  - Lafin  . 


13  n 


n=l 


n 


(16) 


These  are  the  inertia  coefficients  of  the  system  (cf.  Ref.  2, 
p.  188). 


III.  Application  to  a  Floating  Rectangle 


If  the  curve  C  is  symmetric  about  the  x  axis,  then 
in  the  conformal  mapping  this  axis  may  be  mapped  into  the  real 
axis  of  the  ^  plane.  In  such  a  case,  a  normal- velocity 
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distribution  which  is  an  odd  function  of  0,  that  is,  one 
for  which 

a  =0 
n 

for  all  n,  corresponds  to  the  condition 

0  =  0  on  y  =  0  . 

This,  however,  is  just  the  condition  for  the  surface  y  =  0 
to  be  a  free  surface.  If,  consequently,  we  are  dealing  with 
a  cylinder  floating  on  the  surface  of  the  liquid  occupying 
the  half-space  y  ^  0  ,  the  projection  of  the  inunersed  por¬ 
tion  of  the  cylinder  being  C^,  then  the  problem  is  equiv¬ 
alent  to  the  previously  treated  one  of  the  submerged  closed 
cylinder,  provided  C  is  the  union  of  and  of  its  reflec¬ 

tion,  and  the  normal- velocity  distribution  is  antisymmetric 
about  the  x  axis. 

The  only  possible  modes  of  motion  are,  therefore,  those 

for  which  the  coefficients  a  vanish.  Furthermore,  since 

n 

the  kinetic  energy  of  the  half-space  is  half  of  that  of  the 
equivalent  full  space,  the  inertia  coefficients  for  the  float¬ 
ing  case  are  half  of  those  given  by  Equation  (16). 

we  are  interested  here  in  the  motion  of  a  rectangular  box 
of  beam  2a  and  draft  b.  The  equivalent  problem  is  that  of  the 
closed  rectangle  given  by 

-  a  <  X  ^  a  , 

-  b  1  y  <  b  • 
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The  six  degrees  of  freedom  which  are  described  in  the  Introduc¬ 
tion  and  illustrated  in  Figure  1  are  conveniently  grouped  into, 

the  three  which  are  antisymmetric,  and  secondly,  the  three 
which  are  symmetric  about  the  ^  axis.  Consequently, 


®i,  em 
s . 

1, em+i 


s 


i,  2m 


®i, 2m+i 


1 

=  4  /  f.(o)sin  2m0(o)da,  m  =  1,2... 
o 


=0,  m=  0,1,2... 


*  0  ,  m  =  1,2. . . 


=  4  /  f . (a)sin(2m+l)0(o)do,  m  =  0,1... 
o  ^ 


i=l,2,3 


i»4,5,6 


It  remains,  then,  only  to  obtain  the  mapping  function  0(o). 
The  appropriate  mapping  is  given  by  the  following  special  case 
of  the  Schwarz-Christof fel  transformation: 

1. 

=  R(i  _  2C"^  cos  2a  +  (17) 


where  R  is  real,  and  the  points  C  *  ±  ®xp(±  1®)  correspond 
to  the  corners  of  the  rectangle.  In  accordemce  with  Equation  (9) 
we  have 


ds 

d0 


R 


2  cos  2a  -  2  cos  29 


(18) 


Equation  (l8)  can  be  integrated  with  the  aid  of  elliptic  in¬ 
tegrals.  Using  the  notation  of  Ref.  3,  we  let 
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k  =  sin  a  ,  k'  =  cos  a  . 


The  coordinates  of  the  rectangle  are  related  to  9  by  the 
equations 


X  *  ±  a, 

y  =  2R[E(k,0)  -  k'®p(k,0)], 

(19) 

X  =  2R[E(k',0«)  -  k®p(k',0')], 
y  =  ±  b, 

\tfhere 

<t>  ■  sin  ^(sin  0/k),  0'  »  sin”^(cos  6/k')  . 

In  particular 4  the  dimensions  of  the  rectangle  are  related  to 
R  and  a  by 

a  -  2Rk'®B', 
b  -  2Rk®B  , 

so  that  the  beam-draft  ratio  2a/b  is  a  function  of  a  only, 
as  is  the  scale  factor  y,  defined  by 


Consequently,  both  a  2md  y  are  functions  of  2a/b.  For 
the  two  ratios  of  2  and  3,  the  values  of  a  and  y,  obtained 
by  interpolation  fr<xn  the  tables  of  Ref.  4,  are  given  in 
Table  1. 
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Since  for  the  closed  rectangle 

p  =  4(a+  b)  , 

the  relation  between  a  and  6  is  given  by 

i 

^  I  2  cos  2a  -  2  cos  29  \  ^  ,  (20) 

Though  Equation  (20)  is,  of  course,  integrable  in  terms  of 
elliptic  integrals,  there  remains  the  task  of  inverse  inter- 
poloation  in  order  to  obtain  6  as  a  function  of  o.  If  this 
task  is  to  be  performed  by  a  digital  computer,  then  it  is  far 
simpler  to  include  the  integration  of  Equation  (20)  in  the 
same  program,  as  well  as  the  calculation  of  cos  n9  and  sin  nd. 
The  results  of  these  computations  are  shown  in  Tables  2  to  3* 

The  inertia  coefficients  for  the  six  aforementioned  degrees 
of  freedom  were  calculated  by  numerical  integration  and  sum¬ 
mation  up  to  n  «  10,  and  are  tabulated  in  matrix  form  in 
Table  6. 
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APPENDIX 


For  the  two  translatory  degrees  of  freedom^  that  is  i»l 
and  i=4  (see  Figure  1),  the  integrals  for  s^^  are  expres¬ 
sible  in  closed  form. 

Consider,  first,  i^l.  We  have 

y  a  1 

s  *  ”  o  /  2mS(2  cos  26-2  cos  2a)  ^  d6  (Al) 

but  the  sine  of  an  even  multiple  of  6  can  be  written  as 


.in  2«e  .  2  .in  0  CO.  9  (^li) 

The  change  of  variable 


sin  6  >  sin  a  sin  ^ 


yields 


(A3) 


Similarly,  for  1«4, 


n/2 


Atzm*x 


j  8in(2ro-i-l)e  (2  cos  aci  -  2  cos  26)*de,  (a4) 


but 


sin(2m  +  l)e  -  (-l)®sin  [ikj^[5-ki[  cos«^6  .  (A5) 
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On  changing  the  variable 


we  have 


cos  0  =  cos  a  sin  ^  ^ 


s 


=  (-ir7 


4>  sm+i 


cos®a 


m 

2 

k»0 


2k„  7t/2 

—  /  cos^V'sin®^ 


m 

cos^a  2 
k-O 


(-l)Nm-Hc)  'cos^^^ 
kl(k+  l)  i(m-l) !  • 


In  the  particular  case  a/b  >  1,  for  which 

co8®a  «  1/2  , 

we  have 


(a6) 


A 

Try 


1, 


n»  1 


2(-l) 


2r-2)  I 
r-l)i  * 


n  »4r-l,  r  ■  1,2, ... 


(A7) 


0 


all  other  n  . 


From  the  closed^form  expressions  (A3)  and  (a6),  the  iner¬ 
tia  coefficients  m  and  m  can  be  ccxnputed  to  any  desired 

11  44 

accuracy,  for  we  can  find,  with  the  aid  of  Parseval's  theorem, 
an  upper  bound  to  the  trxincation  error  due  to  summing  a  finite 
number  of  terms  on  the  right-hand  side  of  Equation  (l6).  De¬ 
fining 
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we  have 


(N) 


IHii 


PP®  ^  °in  *  ®in 
~  n-1  ^ 


(N) 
Am.  . 

XI 


(N)  oo®  *  +8 


mii-mii 


P£ 

n 


z 

n»N+l 


in  xn 
n 


Hence 


(N)  OD®  “ 

^®ii  1  7r(N+  1)  1  ^®in  ®in^ 


n-M+1 


(a8) 


but 


(A9) 


Calling  the  integral  on  the  right-hand  side  of  Equation  (A9) 
we  have 


('In*  •!«)!• 

n“i 


(10) 


In  particular,  for  the  cases  under  consideration 

a 

K  -  4(7/8)=/  (2  cos  26-2  cos  2o)  d6 


and 


(7/4)= (sin  2a  -  2a  cos  2a), 


n/2 

K  -  4(y/8)=  /  (2  cos  2a  -  2  cos  26)  d6 

*  a 

■  (7/4)=[sin  2a  +  2(J  -  a)  cos  2a). 


-15- 


TABLE  1 


Mapping  Parameters 


2a/b 

a 

y 

2 

0.785398 

1.180341 

(-ttA) 

3 

0.693231 

1.175372 
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Washington  25,  0.  C. 

Attni  Ad3 

Rs 

Ra3 

Rs3 

RsN 


'1; 

,4  ' 

1/ 

1 

,1  I 


Director  of  Naval  Intelligence 
Navy  Departammt 
Washington  25,  D.  C. 

Attni  OP-922V  (1) 


Chief,  Bureau  of  Yards  and  Dodcs 
Department  of  the  Mmvy 
Washington  25,  D.  C. 


Attni  Cods  D-213  (1^ 

Code  D-222  ,1, 

Code  D-410C  ,1 

Code  D-440  ,1 

Coda  D-500  (1 


MAVY  (Continued) 

Commanding  Officer  and  Director 


David  Taylor  Hodel  Baain 
Washington  J,  D.  C. 

Attni  Code  UO  (1) 

Code  600  (1) 

Code  700  (1) 

Code  720  (1) 

Code  725  (1) 

Code  731  (1) 

Code  740  (2) 

Commander 

U.  8.  Naval  Ordnance  Laboratory 
White  Oak,  Maryland 
Attni  Tech.  Library  (2) 

Tech.  Evaluation  Dep.  (1) 
BE  (1) 

EH  (1) 

R  (1) 

Director 

Material  Laboratory 
Haw  York  Maval  Shipyard 
Brooklyn  1,  Hew  York  (1) 

Commanding  Officer  and  Director 
U.  8.  Haval  Elactronica  Lab. 

San  Diego  52,  California 
Attni  Code  4223  (1) 


Officer- in-Charge 
Haval  Civil  Bngin.  Research 
Haval  Civil  Bngin,  Research  and 
Evaluation  Laboratory 
U.  8.  Haval  Construction 
Battalion  Canter 
tort  Husneme,  California 
Attni  Coda  753  (l) 

Director 

Haval  Air  Experimental  Station 
Haval  Air  Material  Center 
Haval  Base 
Hiiladelphia  12,  Pa 
Attni  Materials  Laboratory  (1) 

Structures  Lab.  (l) 

Officer- in-Charge 

Dkiderwatar  Ej^losion  Research  Div. 

Horfelk  Haval  Shipyard 

Portasnuth,  Virginai 

Attni  Or.  A.H.  Kail  (2) 


Commander 

U.S.  Haval  Providing  Grounds 
Dahlgren,  Virginia  (1) 

Super intendent 
Haval  Gun  Factory 

Washington  25,  0.  C.  (1) 

Commander 

Haval  Ordnance  Test  Station 
Inyoksm,  China  Lake,  California 
Attni  Physics  Division  (1) 

Mechanics  Branch  (l) 

Commander 

Haval  Ordnance  Test  Station 

Underwater  Ordnance  Division 

3202  E.  Foothill  Boulevard 

Pasadena  8,  California 

Attn:  Structures  Division  (l) 

Commanding  Officer  and  Director 
Haval  Engineering  Experiment  Station 
Annapolis,  Maryland  (1) 

Supe  r in  tendent 

Maval  Post  Graduate  School 

Monterey,  California  (1) 

Commandant 

Marine  Corps  Schools 

Quantico,  Virginia 

Attni  Director,  Marine  Corps 

Development  Center  (1) 

AIR  gORCB 


Commanding  General 
U.  S.  Air  Force 
Washington  25>  D.  C. 

Attni  Research  and  Developawnt  Div. (1) 


Commander 

Air  Materiel  Command 
Wright-Patterson  Air  Force  Base 
Dayton,  Ohio 

Attni  WOOSZ  (1) 


Commander 

U.  S.  Aj.r  Force  Institute  of  Technology 
wright-Patterson  Air  Force  Base 
Dayton,  Ohio 

Attni  Chief,  Applied  Nech.  Group  (l) 


Director  of  Intelligence 
Headquarters,  D.  S.  Air  Force 
Washington  25,  D.  C. 

Attni  P.  V.  Br.  (Air  Targets  Div.)  (1) 
AIOZH-ia2  (2) 


Commander 

Air  Reaearch  and  Development 
Command 
P.  O.  Box  1395 

Attn:  RDNPE  (l) 

Commander 

NBDD 

Wtight-Patteraon  Air  Force  Baee 
Oaytonj  Ohio 

Attn:  WNRC  (1) 

miMIDS  (1) 

NNRMDD  (1) 

Commanding  Officer 

USMMOBU 

Kirtland  Air  Force  Baee 
Albuquerque,  Hew  Mexico 
Attn:  Code  20 

(Or.  J.  H.  Brennan)  (l) 

0TH8R  OOVBRXBWT  ACTIVITIB8 

U.  8.  Atomic  Energy  Commiaaion 
Maahington  25,  0.  C. 

Attn I  Director  of  Reaearch  (2) 

Director 

National  Bureau  of  Btandarda 
NRahington  25,  0.  C. 

Attn:  Oivialon  of  Nechanica  (l) 
Bngin.  Nechanica  See.  ll) 
Aircraft  Structurea  (l) 

CooBBander 

U.  8.  Coaat  Guard 

1300  1.  Street,  H.  N. 

Naahington  25,  D,  C. 

Attn I  Chief,  Teating  and 

Development  Div.  (1) 

U.8.  Maritime  Adminiatration 
General  Adminiatration  Office  Bldg. 
Maahington  25,  O.  C. 

Attn I  Chief,  Div.  of  Preliminary 
Oeaign  (1) 

National  Adviaory  Committee  for 
Aeronautica 
1512  B.  Street,  N.M. 

Maahington  25r  C. 

Attn:  Loada  and  Structurea  Div. 

(2) 

Director 

Langley  Aeronautical  Laboratory 

Langley  Field,  Virginia 

Attn:  Structurea  Div.  (2) 


Director 

Foreat  Producta  Laboratory 

Madiaon,  Miaconain  (l) 

Civil  Aeronautical  Adminiatration 
Department  of  Commerce 
Maahington  25,  D.  C. 

Attn:  Chief,  Airframea  and  Equip. Br.(l) 

National  Sciencea  Foundation 
1520  H.  Street,  N.  M. 

Naahington  25,  D.  C. 

Attn:  Engin.  Science  Diviaion  (l) 

National  Academy  of  Science 
2101  Conatltution  Aventae 
Naahington  25,  D.  C. 

Attn:  Tech.  Director,  Committee  On 

Shipa'  Structural  Design  (l) 

Executive  Secretary,  Conmittaa 
on  Underwater  Mar  fare  (l) 

Director,  Operations  Reaearch  Office 

John  Hopkina  University 

7100  Connecticut  Avenue 

Chevy  Chaae,  Maryland 

Naahington  I5,  D.  C.  (l) 

Or.  Alvin  C.  Graves,  Director 
J-Division,  Loa  Alamos  Scientific  Lab, 

P.  O.  Box  1663 

Loa  Alamoa,  Mow  Mexico  (l) 

U.  8.  Atomic  Energy  Commision 
Classified  Technical  Library 
1901  Constitution  Avenue,  H.  M. 
Maahington,  D.  C. 

Attn:  Hra.  Jean  M.  O'Leary 


for  Or.  Paul  C.  Fine  (l) 

U.S.  Atomic  Energy  Commission 
Claaaifled  Tech.  Library 
Tech.  Information  Service 
1901  Constitution  Avenue,  N.M. 

Attni  Mrs.  Jean  N. O'Leary  (1) 

Sandia  Corporation,  Sandia  Base 

Albuquerque,  New  Mexico 

Attn:  Or.  Nalter  A.  NacNair  (1) 

Legislative  Reference  Service 
Library  of  Congress 
Maahington  25,  0.  C. 

Attn:  Dr.  B.  Nahk  (l) 

U.S.  Atomic  Energy  Commission 
1901  Constitution  Avenue,  N.M. 
Naahington,  D.  C.  (1) 


OTHER  QOVERMMEOT  ACTIVITIES  (Contlnuad) 

Massachusatts  Instltuta  of 
Tochnology 

Cambrldgo  39,  Masaachuaatts 
Attni  Or.  Char las  H.  Morris  (l) 

INVEST lOATORS  ACTIVELY  ENOAGEO 

ff  WWTIP  BB8BMWP _ 

Profaator  Lynn  s.  Baadla 
Fritz  Englnaarlng  Laboratory 
Lahlgh  Unlvaralty 
Bsthlaham,  Pannsylvanla  (1) 

Profaasor  H.  H.  Blalch 
Dapartaant  of  Civil  Englnaarlng 
Coluabla  Unlvaralty 
Haw  York  27,  Maw  York  (1) 

Profaasor  R.  L.  Blspllnghoff 
Dapt.  of  Aaronautleal  Englnaarlng 
Nassacihuaatts  Instltuta  of  Tadmo. 
Caabrldga  39,  Massaehusatts  (l) 

Profaasor  B.  A.  Bolay 
Oapartawnt  of  Civil  Englnaarlng 
Coluabla  unlvaralty 
Maw  York  27,  M.Y. 

Profaasor  Bugana  J.  Brunalla,  Jr. 
OapartaMnt  of  Aaronautleal  Engln. 
Prlneaton  University 
Prlncaton,  Maw  Jarsay  (1) 

Profaasor  O.F.  Carrlar 

Flares  Hall  * 

Harvard  University 

Caabrldga  38,  Massaehusatts  (1) 

Profassor  J.B.  carmak 
OepartsMnt  of  Civil  Engln. 

Colorado  state  University 

Fort  Collins,  Colorado  (1) 

Profassor  Herbert  Osraslawlez 
Dapartaant  of  Civil  Engln. 

Coluaibia  University 

632  Mast  125th  street 

Hew  York  27,  M.Y.  (1) 

Profassor  Lloyd  Donnell 
Dapartatant  of  Mechanics 
Illinois  instltuta  of  Taduiology 
Tadhnology  Canter  ^  ^ 

Chicago  I6,  Illinois  (1) 


Profassor  D.  C.  Druckar,  Chairman 
Division  of  Engineering 
Brown  University 

Providanca  12,  Rhode  Island  (1) 

Profassor  A.  C.  Erlngan 

Dapartment  of  Aeronautical  Engineering 

Purdue  University 

Lafayette,  Indiana  (l) 

Profassor  w.  Flugga 

Department  of  Mechanical  Englnaarlng 

Stanford  University 

Stanford,  California  (1) 

Nr.  Martin  Ooland 

Midwest  Rasaard)  Institute 

4049  Pannsylvanla  Avenue 

Kansas  City  2,  Missouri  (1) 

Profassor  J.  M.  Qoodlar 
Department  of  Naehanleal  Engineering 
Stanford  University 

Stanford,  California  (1) 

Profassor  L.E.  Goodman 
Bnglnaar  Bx^rlmant  Station 
University  of  Minnesota 
Minneapolis,  Minnesota  (1) 

Profassor  W.J.  Hall 

Dapartaant  of  Civil  Bnginasring 

University  of  Illinois 

Urbana,  Illinois  (1) 

Profassor  R.P.  Barrington,  Hoad 
Department  of  Aeronautical 
Engineering 

University  of  Cincinnati 
Cincinnati  21,  Ohio  (1) 

Profassor  M.  Hotanyl 

The  Technological  Institute 

Morthwestarn  University 

Evanston,  Illinois  (1) 

Professor  P.G.  Hodge 
Department  of  Nechanies 
Illinois  Institute  of  Technology 
Technology  Center 

Chicago  16,  Illinoia  (1) 

Profassor  H.J.  Hoff 
Stanford  university 
Stanford,  California 


(1) 


INVESTIGATORS  ACTIVELY  ENGAGED 
IN  RELATED  RESEARCH  (Continued) 
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« 


Professor  N.  H.  Koppmann,  II 
Department  of  Mechanical  Engr. 
John  Hopkins  University 
Baltimore,  Maryland  (1) 

Professor  Bruce  G.  Johnson 

University  of  Michigan 

Ann  Arbor,  Michigan  (1) 

Professor  J.  Kempner 
Department  of  Aeronautical  Engr. 
and  Applied  Mechanics 
Polytechnic  Institute  of  Brooklyn 
99  Livingston  Street 
Brooklyn  2,  Mew  York  (1) 

Professor  H.  L.  Langhaar 
Department  of  Theoretical 
and  Applied  Mechanics 
University  of  Illinois 
Urbana,  Illinois  (1) 

Professor  B.  J.  Lazan,  Director 
Engineering  Experiment  Station 
University  of  Minnesota 
Minneapolis  14,  Minnesota  (1) 

Professor  E.  H.  Lee 

Division  of  Applied  Mathematics 

Ero%m  University 

Providence  12,  Rhode  Island  (1) 

Professor  George  H.  Lee 
Director  of  Research 
Rensselaer  Polytechnic  Institute 
Troy,  New  York 

Mr.  N.  N.  Lamcoe 

Southwest  Research  Institute 

8500  Culebra  Road 

San  Antonio  6,  Texas  (1) 

Professor  Paul  Lieber 
Geology  Department 
Rensselaer  Polytechnic  Institute 
Troy,  New  York  (1) 

ProfEssor  Hsu  Lo 

School  of  Engineering 

Purdue  University 

Lafayette,  Indiana  (1) 

Professor  R.  0.  Nindlin 
DepartsMint  of  Civil  Engineering 
ColvMd>ia  University 
Mew  York  27,  M.  Y.  (I) 


Dr.  A.  Madai 

136  Cherry  Valley  Road 

Pittsburgh  21,  Pennsylvania  (1) 

Professor  William  A.  Mash 
Department  of  Engineering  Mechanics 
University  of  Florida 

Gainesville,  Florida  (1) 

Professor  N.  M,  Metenark 

Department  of  Civil  Engineering 

University  of  Illinois 

Urbana,  Illinois  (1) 

Professor  Aris  Phillisps 
Oepartstent  of  Civil  Engineering 
15  Prospect  Street 
Yale  University 

Mew  Haven,  Connecticut  (1) 

Professor  w.  Prager,  Chainaan 
Physical  Sciences  Council 
Brown  University 

Providence  12,  Rhode  Island  (1) 

Professor  E.  Reissner 
Department  of  Mathematics 
Massachusetts  Institute  of  Technology 
Cambridge  39,  Massachusetts  (1) 

Professor  M.  A.  Sadowsky 
Department  of  Mechanics 
Rensselaer  Polytechnic  Institute 
Troy,  How  York  (1) 

Dr.  B.  N.  Shaffer 

Departswnt  of  Mechanical  Engineering 
How  York  University 
tS  Fourth  Avenue 

How  York  53.  How  York  (1) 

Professor  C.  B.  Smith 

College  of  Arts  and  Science 

Department  of  Nethomatics 

weiker  Ball 

University  of  Florida 

Gainesville,  Florida  (1) 

Professor  J.  Stallswyer 
Departswnt  of  Civil  Engineering 
University  of  Illinois 
Urbana,  Illinois  (1) 

Professor  Eli  Sternberg 
Brown  university 

Providence  12,  Rhode  Island  (1) 


INVESTIGATORS  ACTIVELY  ENGAGED 
IN  RELATED  RESEARCH  (Continued) 


Professor  S.P.  Timoshenko 
School  of  Engineering 
Stanford  University 
Stanford,  California  (1) 

Professor  A. A.  Velestos 
Department  of  Civil  Engineering 
University  of  Illinois 
Urbana,  Illinois  (l) 

Professor  Enrico  Volterra 
Department  of  Engineering  Hech. 
University  of  Texas 
Austin  12,  Texas  (1) 

Professor  Dana  Young 

Yale  University 

New  Haven,  Connecticut  (l) 

Project  Staff  (10) 


